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Abstract 

Turbulence may have been produced in the early universe during several kind of non-equilibrium 
processes. Periods of cosmic turbulence may have left a detectable relic in the form of stochastic 
backgrounds of gravitational waves. In this paper we derive general expressions for the power 
spectrum of the expected signal. Extending previous works on the subject, we take into account 
the effects of a continuous energy injection power and of magnetic fields. Both effects lead to 
considerable deviations from the Kolmogorov turbulence spectrum. We applied our results to 
determine the spectrum of gravity waves which may have been produced by neutrino inhomogeneous 
diffusion and by a first order phase transition. We show that in both cases the expected signal may 
be in the sensitivity range of LISA. 

1 Introduction 

The high isotropy of the cosmic microwave background radiation (CMBR) testifies that the universe 
was a quiet child at photon decoupling time. This does not prevent it, however, to have been a quite 
turbulent baby! 

Several violent and interesting phenomena may have taken place before the universe became mat- 
ter dominated. Phase transitions, and re-heating at the end of inflation are examples (we will discuss 
another one below) of processes which may have taken the universe through a phase during which 
thermal equilibrium and homogeneity can be temporarily lost. Turbulence is expected to have devel- 
oped under such conditions due to the huge value of the Reynolds number. It would be extremely 
interesting for cosmologists to detect any observable relic of such turbulent periods as this may shed 
light on the early stages of the universe and on fundamental physics which is not yet probed in the 
laboratories. Unfortunately, due to the finite thickness of the last scattering surface, density or metric 
perturbations which may have been produced by any causal process before the matter radiation equal- 
ity time than can hardly give rise to detectable imprints on the CMBR. The detection of cosmological 



background(s) of gravitational waves (GWs) may offer the only possibility to probe turbulence in the 
remote past. Although any direct detection of GWs is still missing, this subject is not any more purely 
academic as a new class of ground and space based observatories of GWs are under construction or 
in advanced project. The amplitude and frequency sensitivity of some of these instruments are in the 
proper range to probe many interesting astrophysical and cosmological processes [||. The most inter- 
esting project from our point of view is LISA (Laser Interferometer Space Antenna) which is scheduled 
to be launched around the end of this decade and will possibly achieve a sensitivity in /igfigw of 10 -12 
at milliHertz frequencies Q. 

So far the interest on GWs produced by cosmological turbulence has been mainly stimulated by the 
possibility that turbulence could have been generated at the end of the electroweak phase transition 
(EWPT) [||, The idea is that, if the EWPT is first order, turbulence should have been produced 

when the expanding walls of bubbles containing the broken symmetry phase, or the shock waves 
preceding them, collide at the end of the transition. As discussed in [Q, ||] the stochastic background 
of GWs produced during this process may be detectable by LISA. 

Recently two of us [|| noticed that another, possibly detectable, background of GWs may have been 
produced before neutrino decoupling, that is much later than the EWPT. This mechanism requires 
that the net lepton number density (JV a (x) = ra„ a (x) — np a (x), a = e,/i,r) of one, or more, neutrino 
species was not uniform and changed in space over some characterist scale which was smaller than 
the Hubble horizon at the neutrino decoupling time. A net flux of neutrinos should then be produced 
along the lepton number gradients when the neutrino mean free path became comparable to the size 
of the fluctuations in the neutrino number. It was shown that, depending on the amplitude of the 
fluctuations, such a flux of neutrinos may be able to stir chaotically the cosmic plasma producing 
magnetic fields and GWs. Fluctuations in the neutrino number, as those which are required to power 
all this rich set of effects, can be a by-product of leptogenesis at the GUT scale (this is possible for 
example in the Affleck-Dine 0| scenario of leptogenesis) or be a consequence of active-sterile resonant 
neutrino conversion before neutrino decoupling ||]. 

The aim of this paper is to estimate the characterists amplitude and spectral distribution of GWs 
produced by the mechanism discussed in Q and during a first order phase transition. Our treatment 
follows that of Kosowsky, Mack and Kahniashvili [Q] (KMK) and extend it by introducing some 
physically relevant generalizations. 

In Sec. ^ we will show how the Kolmogorov spectrum of turbulence (which is that adopted by 
KMK) is modified when energy injection into the primordial plasma takes place over a continuos 
range of scales. This is motivated by the observation that neutrino number fluctuations, which might 
power turbulence, generally will not have a unique size but have a more complex spectral distribution. 
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The same consideration can be naturally applied to other kind of stochastic sources. Since strong 
magnetic fields are expected to be produced by the mechanism presented in Ref. |6| and during first 
order phase transitions j|, in Sec. || we will also investigate possible Magneto-Hydro-Dynamical 
(MHD) effects on the turbulent spectrum. We will see that these kind of effects may give rise to quite 
drastic consequences on the expected GW signal. In Sec. || we apply the results of the previous two 
sections to determine the general expression of the GW power spectrum produced by cosmological 
turbulence. In Sec. ||| we estimate the signal to be expected at LISA because of GWs generated by 
neutrino inhomogenoeus diffusion. In Sec. ^ we will determine the characteristics of the GW signal 
produced by a first order phase transition paying particular attention to the case of the EWPT. We 
will show that MHD effects may play a crucial role by enhancing the expected signal. Finally Sec. [7| 
contains our conclusions. 



2 Turbulence spectrum from stirring at any scale 

Three-dimensional turbulence shows a cascade of energy from larger length scales to smaller ones 
(direct cascade): largest eddies form at the scale at which the stirring source acts and after a few 
revolutions they break down into smaller eddies; this process goes on until the damping scale is reached 
(see below). Call £k the rate at which a (momentum) scale k receives turbulent energy (normalized 
to the enthalpy w = p + p of the fluid) 

1 (dpk)m (dul)- m 
£k ~(p + p) dt ~ dt ' lj 

where Uk is the typical fluid velocity at the scale k, and we assumed constant energy density p and 

pressure p, i.e. that the fluid is incompressible. If the process is (almost) stationary, must equal 

the rate at which energy is transferred to smaller scales, i.e. 

St ~ u\k , (2) 

where we assumed (as it is usually done) that the rate at which an eddy breaks down into smaller 
ones is roughly its turnover angular frequency lo^ ~ ku^. 

In general stirring of the fluid by an external source does not take place at a single scale. Rather, 
turbulence in the early universe is likely to be produced by stochastic sources injecting energy over 
a finite range of wavenumbers. Consider any momentum k within this interval. An infinitesimally 
larger momentum k + dk, apart from the energy from the turbulent cascade, will receive an additional 
external power P ext (k)dk, due to the external stirring at k + dk, that is 

"it - F «< fc » ■ < 3 > 
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This is different from conventional Kolmogorov type turbulence where energy is supposed to be injected 
at a single scale and e k = e is momentum independent. 

Following a standard approach (see e.g. [l^] for a pedagogical introduction to turbulence) we 
define the energy spectrum E{k) by 

POO 

(u 2 (x)) EE / E(k)dk , (4) 
J 



with (• • •) standing for ensemble averages; by combining eqs. (0-[|), and noticing that eddies satisfy 
the "uncertainty" relation Ak/k ~ 1 [], we find the spectrum 



E(k) 



k P ext (k')dk' 

k s 



-| 2/3 

/r 5/3 , (5) 



where fcs is the largest length scale at which we inject energy in turbulent motions. Note that eq. (|B|) 

reduces to the usual Kolmogorov spectrum for a delta-like stirring spectrum P ex f 

The energy cascade stops at the damping scale kp, characterized by a local Reynolds number 

Re/) ee U£> of order unity, i.e. 
k D rj 

kD ~ e^t T 3/4 , (6) 

where £tot = / Pext{k)dk is the total external stirring power. The energy spectrum E{k) has ku as 
"ultraviolet" cutoff since any energy injected at k > krj is dissipated into heat. 

We expect substantial deviations from Kolmogorov spectrum if the integral appearing in eq. (||) 
is substantially increasing with k; in such cases we parametrize it as a power law function of k for 

k > k s , 

e k = / P ext (k')dk' = s s ( — 
Jk s \ ks 

with 7 > 0; this gives | 



37/2 

(7) 



E(k) ~ 4 /3 kg 1 /t 7 " 5/3 . (8) 
Again, Kolmogorov spectrum would correspond to 7 = 0. 

It is helpful to rewrite the spectrum and all the other quantities parametrizing them with the 
largest-scale characteristic velocity us rather than with the characteristic rate £5. This permits a 
more intuitive and direct generalization to the relativistic case (see below). From eqs. (Q) and ([/]) we 
obtain 

'£A 1/3 = u (ky /2 - lpi 



u k ~[ T ) = u s[ v J , 0) 



1 Properties of Fourier transforms automatically include the uncertainty relation AxAk ~ 2-k; moreover, it is clear 

that the spatial width Ax of an eddy is roughly its characteristic wavelength A = 2-k jk. This gives Ak/k ~ 1. 

2 Eqs. ([]) and (^) have to be intended as a rough estimate of the real physical spectrum. Although they should give 

a good approximation for small values of 7 (which corresponds to some sort of adiabaticity in momentum space) they 

are probably very crude for 7 of the order of unity. We thank G. Falkovich for making us this comment. 
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where us = (es/^s) We thus have 



U k ~ U S 4 /3_7/2 ^ 2+2 / 3 (10) 

s(jfe) ~ 4^ /3 - 7 f- 5 / 3 (ii) 



1/3-7/2 

fc */o-7/. ^ ^5 _ (12) 



- 4/3-7/2 us k s 



V 

Notice that for 7 > 2/3 the characteristic velocity u^ increases with k: this means that even if 
on large scales the plasma has non-relativistic turbulent velocities, it will reach relativistic velocities 
at a certain kji, if the latter is smaller that the damping wave number kp. In such a case the 
characteristic turbulent velocity for k > kn will be saturated at ~ 1 Q clearly this will be the 
most interesting case for the production of gravitational waves. In order to extend to that relativistic 
case all the formulae we will derive throughout the paper, it is sufficient to replace 7 with 2/3 and us 
with 1: this is straightforward to understand by noticing that with those formal replacements eq. (||) 
produces a costant velocity spectrum with u^ equal to 1. In particular, by applying that prescription 
to Eqs. tfnHiJ) we obtain 

{u k ~ k 
E{k) ~ k- 1 ■ (13) 
k D ~ 1/?? 

A final remark to this section is in order. The standard theory of turbulence, on which we based 
our previous considerations, has been formulated and tested only for non-relativistic fluid velocities. 
A relativistic theory of turbulence is, unfortunately, not yet available. Therefore, it is fair to say that 
by applying our results to situations with relativistic velocities we do a somewhat blind extrapolation. 
We think, however, that our results provide a correct order of magnitude estimate of the real physical 
quantities. The same attitude was adopted by the authors of previous works on the subject P, 0]. 



3 Velocity and magnetic field spectrum for MHD turbulence 

Turbulence in an electrical conducting fluid in the presence of magnetic fields is known as MHD 

turbulence. The dynamics of MHD turbulent energy cascade may differ considerably from that of 

conventional, Navier-Stokes (NS) turbulence, which we discussed in the previous section. The main 

3 More precisely, Uk will saturate to the sound velocity c s which is 1/V3 for a relativistic fluid. The excess of energy 
injected into the fluid will be converted into supersonic shock waves. Although the approximation to treat the fluid as 
incompressible clearly fails in this case, it is reasonable to assume that the presence of shock waves will only make more 
involved the energy injection process without affecting turbulence properties. Notice that ut will be always smaller than 
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reason of such different behaviour resides in the so called Alfven effect |12j : in the presence of a strong 
background magnetic fields Bq the velocity and magnetic field fluctuations become strongly correlated, 



i.e. v ~ ±5B/^/47r(/9 + p). Fluctuations of this kind correspond to linearly polarized Alfven waves 
propagating along the field Bo with velocity 

the so called Alfven velocity. In the absence of external magnetic fields Bo has to be intented as the 
average field computed over a scale Lq L = 2ir/k . In our case we may assume Lq to coincide 
with Ls = 2n/ks- MHD corrections to the turbulent cascade can only be neglected if the turbulent 
velocity u k is larger than va- In the opposite limit MHD turbulence sets-in which may be figured as 
an ensemble of stochastic Alfven waves. 

Let's define k eq as the wavenumber at which u k first becomes equal to va- In the case u k is a 
decreasing function of k, i.e. 7 < 2/3, we may assume a NS type spectrum for k < k cq and determine 



k eq by equating eq. (H) to (14). It is evident that MHD turbulence will never develop if k cq > ko, 



where kn is the damping momentum for the NS turbulence. In the opposite case, however, there will 
be a range of momenta ko > k > k eq where we may expect noticeable corrections to the expressions 
written is Sec. ^. Indeed, since the Alfven velocity is larger than vJ^ s in this range, the time of 
interaction of two eddies with similar momentum Ta ~ 1/kvA is smaller than the turnover time T k . 
Since many interaction events are necessary to change the wave packet amplitude appreciably, the 



energy transfer time becomes longer [12 



T k - T' k = T A ( = 3- » T k ■ (15) 



,TaJ T a 

This effect leads to a milder turbulent spectrum. In the case of stationary energy injection localized 
to a single momentum ks, eq. ( |15| ) implies 

e k ~ u\— (16) 

VA 

E{k) ~ (eva) 1 / 2 k- 3 / 2 (17) 

which is the so called Iroshnikov-Kraichnan spectrum. In the case of a continuous external stirring 
power P ext and 7 < 2/3, ensuring a decreasing slope of the velocity u k , eq. (||) generalizes to 

E(k) ~ (£ S v A ) 1/2 kp y kUi- 1 ) . (18) 

Since the magnetic fields at the GW emission time is not a direct observable, it is convenient to express 
the previous results in terms of k eq which may be directly observed in the GW power spectrum. 
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From our previous considerations it follows that k eq is defined by 

u\k cq ) = e 2 J 3 kpk£* = v\ . (19) 



By substituting this expression in eq. ( |i8[ ) we get the result 



2 / 3 £.-71,7-5/3 k <r k 

g rug ru j ru \ A/gq 



S o kg kec[ ® k ^ \ 2 ) , k ^> k C q . 

Since the GW power spectrum is directly related to E(k) (see section ||), we can already anticipate 
here that possible observations of a variation in the the slope of the GW power spectrum, fitting the 
power laws reported in eq. ([H]), would provide information about the very presence and the strength 
of a smooth (non-turbulent) magnetic field component at the GW emission time. 

Since Alfven waves carry the same amount of magnetic and kinetic energy, approximate equipar- 
tition between kinetic and magnetic turbulent energies should hold for k > k eq . According to the 
authors of Ref. fl3|| , tangled magnetic fields will also contribute to GW emission. Furthermore mag- 
netic fields may have produced observable effects on the CMB anisotropies jL3|, [jl]] and provided the 
seeds of galactic and intergalactic magnetic fields ju]]. Therefore it is worthwhile to write here also 
the magnetic field turbulent fluctuation spectrum. By neglecting the turbulent magnetic energy above 
the equipartition scale we have 

!0 « k A^pn 

2/3.-7 J-* (21) 
Eg kg koQ ^ /c 2 V 2 ) ^ k ^> A^gq , 

where EM{k) is the spectrum of energy stored in magnetic field fluctuations per unit enthalpy of the 
fluid. 

In order to cast the spectrum ( |20| ) in a more symmetric form, we define a new parameter 5 and a 
new scale k' s 

5 = \l + \, k' s - 5 = k^ki~^, (22) 

E{k)~\ flYJ^ k<keq (23) 



so that eq. (|20|) becomes 



The same procedure with the same 5 and k' s works for the relation between a turbulent scale k and 
the corresponding characteristic frequency because ~ ^k 3 ^ 2 E(k) 1 / 2 , see eqs. (1C) and ([TT|). 
Thus we have 

( 1 1/3^-7/2^2/3+7/2 k<k 
f, ~ ) 2tt 6 6 eq , . 

Jk ~ I JL .I/ 3 y - 5/2 k 2 / 3 + 5 / 2 k>k 

I 2tt £ S K S K > K > K eq . 
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In section |I| we will compute the amount of relic GWs from NS turbulence. The only specific features 
of NS turbulence that we will use are the turbulent spectrum E(k) and the relation between frequency 
and length scale. With the above results in mind, this means that all the formulae of section [| are 
straightforwardly applicable to MHD turbulence by the formal replacements 7 —* 5 and kg — > k' s . 

For 7 > 2/3 MHD corrections will be relevant only if us < va 0- It is easy to verify that all the 
previous expressions which were derived for 7 < 2/3 can be applied to the 7 > 2/3 case by simply 
inverting their range of validity with respect to k. Namely, we get a MHD turbulent spectrum for 
k < k cq and, if k cq < ko, a NS spectrum for k > k eq - 



4 Relic gravitational waves 

We now want to compute the amount of GWs produced by cosmological turbulence. We are interested 
in the present-day energy spectrum normalized to the critical density, 

s 7M ■ (25) 

We will partly make use of the recent results of Kosowsky et al. these were derived under the 
assumptions that the turbulence has a Kolmogorov spectrum (7 = 0, in the notations of section ||) 
and that it lasts for a time interval r which is very short compared to the Hubble time at that epoch. 
We will perform a more general calculation, including the possibility that the plasma is stirred at any 
scale and develops a general turbulent spectrum, as discussed in section We also consider the case 
that the time interval during which turbulence is active is long and comparable to the Hubble time. 
We deal with this case assuming that the coherence time of turbulence, which will be of the order of 
the characteristic turnover time of a single turbulent eddy, is substantially smaller than the Hubble 
time, so that all the machinery developed in the small r approximation is applicable, because different 
eddies act as uncoherent sources of GWs. The latter requirement is automatic if we consider eddies 
whose characteristic length scale is well inside the horizon and whose characteristic velocity is not too 
small. 

Our convention for the Fourier transform of any function F(x) at fixed time is 

F(k) = i Jd 3 x e +ikx F(x) , (26) 

where x and k are physical quantities and we keep the reference volume V for dimensional convenience. 

Let us consider a turbulent fluid in the Friedmann-Robertson- Walker (FRW) Universe at a fixed 
time t* and denote its velocity field as u(x). The turbulent motions of the fluid will produce a 
4 Clearly this will never be the case if us is close to the speed of light. 
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stochastic background of GWs. We define the perturbed metric as g^ v = g^yit) + h^(x,t) where 
g^ v is the FRW unperturbed metric. Under the assumptions that the fluid is incompressible and 
that turbulence is statistically homogeneous and isotropic, the velocity correlator of the fluid can be 
parameterized as 

(«i(k,t) u*(k',t) ) = ^-{Sij - kk) ^ ( k " k ') p ( k )- (27) 
J V 

It is straightforward to show that the velocity spectrum P(k) is related to the turbulent energy 
spectrum E(k) defined in section ^ by 

E (k) = _J_%£^ = }Lk 2 P(k) . (28) 
p + p dk 7r 2 w y J 

If the source acts during a short time interval At* (see the discussion above), much smaller than the 
Hubble time at £*, we can neglect the expansion of the Universe, and use the result of Kosowsky et 
al. for the correlator of GWs produced during the interval At*, 

( MM*) K 3 (k',Q > ~ VV2(l6nGfAUw* §3{k _ y ^ p((?)p(|k _ q|) ? (2g) 

where u? = p + p is the enthalpy density of the fluid 0. 

For a power-law spectrum P(A;) = A k~ a in the range ks < k < kz>, the integral appearing in 
eq. (|2^) is well approximated by M 



1 



d 3 qP(q)P{\k-<i\) ~4ttA 2 



k 3 ~ 2a a + k 3 p 2a k- a k 3 s ~ a 



(30) 



(3-a)(3-2a) 3 - 2a (3 - a) 

for ks < k < ki) and outside that interval. Eq. ( ^0[) holds for a > 3/2: for a < 3/2 the integral must 
be computed more carefully, but this case will not be of interest for us. There are three qualitatively 
different regimes of eq. (30): for a > 3 the last term of the r.h.s. is dominant, while for a < 3 the first 
dominates. The limiting case a — ► 3 gives a term which behaves like A; -3 log k, coming from the sum 
of the first and the last term. Taking into account the relations a = 11/3 — 7 an d AV/ir 2 = u 2 s k 2 J 3 7 , 
coming from comparing eqs. ( |28| ) and ([ll]), we obtain 

f ^ kf'^ , 7 < 2/3 

/ A P(5)P(|k - q|) * ^ 4 fcf^ x (37 3 ^;-f_) 67) fc-»AH*r ? 7 > 2/3 (31) 

I fc- 3 logA, 7 ^ 2 /3 



5 The r.h.s. of eq. (^9|) is 4 times larger than the corresponding quantity appearing in ref. Q. This is because our 
definition of the metric perturbation hij is twice larger than the one used in that paper, where hij = \8gij. Although 
the two descriptions are equivalent, one must pay attention when applying them to calculate the energy density of the 
GW background. In ref. GJ this is done by means of formulae taken from ref. which instead uses our definition of 
hij. Therefore we think that the authors of ref. M] underestimates the GW energy density by a factor of 4. 
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as leading contributions. 

It is now straightforward to compute the energy spectrum of gravitational waves /ig^gw(/) for the 
three different cases of eq. (|3l|). In the Appendix we present all the details of the computation. Here 
we summarize only the strategy and the results. From eq. ( p9| ) one computes the real-space correlation 
function 



( hij(?£, i*) /ijj (x, i*) ) 



V- 



and from that the characteristic amplitude h c (f) of GWs at frequency /, defined by (see ref. 
details) 



( hij(-x,t) hij(x,t) ) 



(32) 
]] for 

(33) 



where one has to express everything in terms of the frequency / of a turbulent scale k by means of 
the relation 

h = g - ± us 4 /3 " 7/2 ^ 2/3+7/2 (34) 

following from eq. (jlOj) . 

The expansion of the Universe redshifts the frequency and damps the amplitude of GWs: both 
are inversely proportional to the scale factor a. If the subscripts and * stand respectively for the 
present time and the time of production, we have [] 



— ~ 1.7 X 10" 
the frequency we observe today is thus 

/o ~ 1.15 x 10~ 7 mHz 



10 


"10.75" 


1/3 


"MeV" 




. 9* . 







MeV 



9* 



10.75 



1/6 



where we used the fact that the Hubble parameter at time of production is 



H* = 1.66 gl /2 — - ~ 4.5 x 10~ 22 MeV 
M P 



MeV 



9* 



10.75 



1/2 



(35) 



(36) 



(37) 



and that 1 mHz ~ 6.7 x 10~ 25 MeV. From eq. (^) we see that in order to have detectable (i.e. high 
enough: for instance LISA operates at / ~ mHz) frequencies we must have turbulent motions at 
length (time) scales well below the horizon and at large temperatures T* > 1 MeV. For such a reason 
we will be interested in the high-frequency region of the spectrum: the maximum frequency we can 
expect today, in terms of the damping scale kr> at time of production, is 



/max = —fk D ^ 2 x 10- 8 mHz I ) 



2/3+7/2 



U S 



ks 1 1/3-7/2 p T> 



MeV 



9* 



nl/6 



10.75 



(38) 



Although we use here a parameterization suitable to describe GW production around neutrino decoupling time, we 
should stress that the expressions contained in this section are completely general. 
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The characteristic amplitude measured today at a frequency / is 



h c (f,to) = —h c (f* = —/,**) , 
ciq a* 



(39) 



and the energy spectrum we are interested in is related to that by p]] 

h c (f) = 1-3 x 1(T 15 v^V(7) • (40) 

Applying this procedure to the three different cases of eq. (|3~l|) leads to (see the Appendix for the 
details) 

2(7-67) 



7 <2/3 

7 >2/3 
7 ~ 2/3 



^0^gw(/) 
^gw(/) 
^0^gw(/) 



^ i n— 4 6 

6 x 10 — —u$ 

^ i n— 4 6 

6 X 10 -j-li^ 



4 x 10~ 5 ^-u? 





-l 


5* 


-1/3 


.-"*. 




.10.75, 






-l 




-1/3 


.-"*. 




.10.75. 






-l 


5* 


-1/3 






.10.75. 





D-, 



f_ 

fs 

f_ 

fs 



4+37 



18(1-7) 
4+3 7 



LY\ og L 

fs) fs 



(41) 
(42) 
(43) 



where C~ 



(11-37) 



(2-3 7 )(4+3 7 ) alld Dr y ~ (37-2) (4+3 7 ) (13-6 7 ) 



are 0(0.1) numerical coefficients and fs is the 



properly redshifted frequency corresponding today to the largest length scale at time of production, 



fs~2x 10- s u s 



ks_ 



MeV 



.9* 



nl/6 



10.75 



mHz 



(44) 



In the Kolmogorov case we get 
7 = =► h 2 n gw (f) ~ 6 x 10" 



32 ^* 


/ mHz N 


\ 7/2 19/2 
) U S 


' k s ~ 


5/2 




7/2 


9* 




V / , 






.MeV. 




.10.75. 



1/4 



(45) 



Apart for the factor 4 which we already discussed, this expression is consistent with the result of KMK 
0| . It is evident from eq. (|4^) that a Kolmogorov turbulent spectrum give rise to an undetectable GW 
intensity at mHz frequencies if T ~ 1 MeV and the energy injection scale is comparable to the Hubble 
horizon size. However, as we will see in more details below, observationally more promising intensities 
can be obtained, even for a Kolmogorov spectrum, for larger emission temperatures especially in those 
case in which the injection momentum ks increases more rapidly than H* with T*. 

For 7 > 2/3 the dependence of h^Qgw on / is milder than in the 7 < 2/3 case: for 7 > 1 the 
spectrum is even increasing with /. For 7 = 1 the spectrum is flat, 
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ZlQ^gW ^ 



10" 



At* 



T a S 



-ff* 



10.75 



1/3 



(46) 
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4.1 Time-integration 



If the time during which turbulence is active is long, i.e. of the order of the Hubble time H~ l , we 
have to perform an integration over production time i* in order to get the correct /iQ^ gw . To be more 
specific, let's consider the case of GWs produced at the time of neutrino decoupling which we will 
illustate in more details in the next section. Consider for simplicity a turbulent spectrum saturated 
at relativistic velocities u^ ~ 1: as discussed in section ^, it is well described by the 7 = 2/3, us = 1 
case. The energy density of GWs produced in a short time interval dt* is given by eq. (|43| ) with the 
us factor set to 1. Neglecting for simplicity the logarithmic factor (this a conservative assumption 



since / > fs) and using eq. (44) we obtain 



h 2 n gw (f) ~ s x 10- 13 f 

Jt, 



8 x 10 



-1.3 



ton 



1 o / mHz \ 
8 < 1(T 13 ( —r- J 



MeV 

10.751 1/6 r T < 



10.751 1 / 6 



9* 



where we assumed T on 3> T q q and g* 



V / 

const. 



9* 
T 

- 1 on 

MeV 



d(T*/MeV) 

^off 

10.751 1/6 



9* 



(47) 



The cutoff frequency is (see eq. (|3q)) 



/max^2x 10- 8 mHz 



MeV 



9* 



10.75 



1/6 



(48) 



For the relativistic case ~ 1 the damping scale is kn ~ 1/r/, see section in the case discussed in 
section [| the viscosity is rj = igrj^x @v — ^^u- The mean free path £ u scales like T~ 5 , and at neutrino 
decoupling it is equal to the horizon length. This gives 



fn 



1 mHz 



40 MeV 



9* 



10.75 



1/6 



(49) 



where we used ~ 1 MeV. 



5 GWs from turbulence produced by neutrino inhomogeneous dif- 
fusion 

In Ref. H two of us showed that if isocurvature fluctuations existed in the early universe under the 
form of a space dependent neutrino net number (N a (x) = n Ua (x) — n Pa (x), a = e, /i,r) turbulence 
should have developed in the primordial plasma before neutrino decoupling. Such an effect arises as a 
consequence of neutrino currents which flow along the lepton number gradients when the neutrino mean 
free path £ U (T) becomes comparable to the characteristic size A of the isocurvature fluctuation. The 
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residual elastic scattering of the diffusing neutrinos onto electrons and positrons will then accelerate 
these particles together with the photons and the rest of the neutrinos to which e~ and e + are still 
tightly coupled (baryons can be disregarded at that time). A random distribution of fluctuations will 
generally give rise to vortical motion of this composite fluid. Depending on the velocity and the size of 
the eddies, turbulence may then develop in the interval of time during which the random forces due to 
neutrino elastic scattering overcome the shear viscosity force [|. It was showed that the fluid velocity 
can approach the speed of light if the amplitude of the fluctuation 5N/N is not too much smaller than 
1. In such a case the authors of || claimed that turbulence should give rise to a stochastic background 
of GWs at a level which may be detectable by LISA space observatory or by its upgrading. In this 
section we investigate this issue in more details. 

The fluid acceleration mechanism is suitably described by kinetic theory in an expanding geometry 



[15|. The first step of the computation is to determine the neutrino momentum flux by solving the 
transport equation for the neutrinos in the presence of a source term given by the fluctuation of N a 
and a collision term due to neutrino-electron(positron) scattering. 

Let's first introduce some useful notation: x = is time in units of the neutrino decoupling time 
td; H(t) = 2^x _1 and Hd = ^ are the Hubble rates for tQCD <t <td and at t = t^; A(i) = 2t ( i\x 1 / 2 
and A = \{tj) Hd is the fluctuation wavelength in natural units and in units of HJ ; k{t) = ^- kx^ 1 / 2 
and k = k(td)/Hd are the corresponding wavenumbers; l v (t) = r = 2td x 5 ^ 2 are the neutrino mean 
free path and collision time; K a (x, i) = i Jk f Ua {E, k) is the specific momentum flux of u a ; 
finally, we define the derivative ( )' = -§-( ) = td4r( )• Then the momentum transport equation can 



Pv a J JV^'^W) 

75S ( ) = td dt 

be written Iql 



! K ^-M^)- 4ffK »->»-^ (50) 

where v is the macroscopic velocity of the e + e~ / y fluid. The Euler equation for v is 

where r\ = jgr^j^ is the shear viscosity due to neutrino diffusion and we assumed for simplicity to 
have the same inhomogeneities in the density of z/s and z^'s, i.e. ~ Kp []. Finally, the evolution of 
the fluctuations 5N a /N a is governed by the continuity equation 

^-N a = -V • (Ar a K a ) . (52) 



7 This is the same of saying that the Reynolds number is much larger than one. 

8 As long as we are not concerned about magnetic field production (see below), we can safely approximate the v — e~ 
and v — e + mean collision time to be equal. The slight differences between the values of l v and r for the and the i^,t, 
are also not significant here. 
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By making explicit the x dependences and passing to the Fourier space we transform the last three 
coupled equations into 

K k ' = -fkS k x- 1 / 2 - 2X- 1 K k - x" 5 / 2 (K k - v k ) 

v' k = 2r x- 5 / 2 (K k - v k ) - ix" 1 v k -%P x 3 / 2 (v k - K k ) (53) 
S' k = -i x -V 2 kK k 

where 5 k is the Fourier transform of (^jr-) and we defined r = . In eq. (|53|) we kept only the 



V N v J p+p 

first order in 5 k and we ignored the Lorentz factor 7 of eq. (^l|) 0. The initial conditions at x = are 
K k = v k = 0, 5 k = 5° k . 

The characteristic decay time x of the fluctuation 5 k can be calculated by exactly solving the 
system (|53|). For x <C x the fluctuation can be considered frozen at its initial value S k ; at that early 
time the system ( |53"| ) reduces to 

K k ' = - i ~k6° k x- 1 ' 2 - 2X" 1 K k - x" 5 / 2 (K k - v k ) 
v' k = 2r x" 5 / 2 (K k -v k )- \x~ x v k , 



where we kept the leading order in 1/x. The system (p4\ ) is easily diagonalized and approximately 
solved by 

K k (x <C x) ~ v k (x <C x) ~ — k5 k x 1/2 , (55) 

6 1 + 5r 

with K k — v k = 0{x 2 ). Plugging this into the continuity equation for 5 k gives x ~ 6 ^^ r ^ as char- 
acteristic decay time. At later times x > x the source term is negligible and the neutrino momen- 
tum K k and the fluid velocity v k are damped by the expansion of the Universe (the viscosity term 
~ k 2 x 3 / 2 (v k — K k ) will become important even later). Therefore the e + e~7 fluid reaches the maxi- 
mum velocity at time x ~ x, 

which has the same k dependence as the initial fluctuation 5 k . This rough analytic estimate is con- 
firmed and made more precise by exact numerical integration of the coupled equations ([53]) : eq. ( [jG| ) 
overestimates the actual value of ^™ ax by a factor 1.2 — 1.3. 

As an illustration, let's first consider GW production in the simple case of isocurvature fluctuations 
peaked at a single comoving wavenumber ks- When l v becomes comparable to ks, let's call this time 
t*, neutrino diffusion give rise to a vortical velocity field according to the mechanism discussed in the 
above. Let's assume for simplicity that the peak velocity v ks = us is reached instantaneously. This 
velocity will be maintained for a time of the order of the Hubble time H^ 1 ~ i*. Since, the velocity 
9 The sound velocity l/\/3 sets an upper bound for macroscopic velocities of the fluid, therefore 1 < 7 < -v/3/2 ~ 1.2. 
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turnover time is of the order kg 1 , which is supposed to be much smaller than H~^, turbulent energy 
cascade from ks down to ko have all the time to develop. This is not enough, however, to conclude 
that fully developed turbulence is actually formed. In order to verify if this is the case we have to 
verify if the Reynolds number is much larger than unity. We have 

Re = ^OtS- ~ 30 u s (57) 
V 

where we used the expression rj = 15 | p p + p ) t v — jjgC for the viscosity and we assumed kg 1 ~ l v . For 
a localized stirring source the standard theory of turbulence would not hold for such a low Reynolds 
number. This is because the required statistical isotropy and homogeneity of the turbulent velocity 
field cannot be achieved in this case. The situation may however be different if the stirring source is 
stochastic itself providing statistical isotropy and homogeneity from the very beginning. We assume 
that in such a situation the results of Secs.0 and can be applied also for order unity Reynolds 



numbers, provided that kp is larger than ks- Furthermore, as follows from eq. (12), kp increases with 
respect to ks when 7 > 0. The most favourable case is that in which isocurvature fluctuations give 
rise to fluid velocity close to unity over a wide range of wavenumbers. As we discussed in Sections |2| 
and ||, this give rise to a saturated turbulent spectrum with 7 = 2/3. In such a case eq. ( fl3"l) holds for 
the GW energy density. 

In general the turbulent velocity power spectrum will depends on the power spectrum of the neu- 
trino number fluctuations which in turn will depend on the mechanism responsible for their generation. 
The computation of such spectrum is beyond the aim of this paper. In Ref. || the authors considered 
two possible scenarios. In the first one N v fluctuations with amplitude of order 1 are produced as a 
consequence of active-sterile neutrino oscillations according to the mechanism proposed by Di Bari 
[|| (see also Ref. [0]). The seeds of the neutrino domains are tiny fluctuations in the baryon num- 
ber which may have been produced during the QCD phase transition or the inflation. Their power 
spectrum will determine that of neutrino domain hence also that of GWs. A continuous spectrum 
is generally to be expected. Although in principle this mechanism may naturally give rise to very 
large amplitude isocurvature fluctuations and GW signal, unfortunately it would be probably hardly 
detectable by LISA. This is a consequence of the low critical temperature at which neutrino domains 
may have been formed which turns into a low GW frequency (see eq. (p9|)). Better observational 
perspectives are offered by the second mechanism considered in || which is based on a generalization 
of the Affleck-Dine baryogenesis mechanism |7], |17[| . In this case neutrino domains may be formed 
during inflation when a scalar field carrying lepton number rolls along (nearly) flat directions of the 
superpotential. Isocurvature fluctuations of amplitude as large as 1 may be formed without invok- 
ing too extreme assumptions. Although the fluctuations power spectrum is expected to be strongly 
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model dependent, a nearly flat spectrum is a quite reasonable possibility. As we discussed above, a 
flat spectrum of the neutrino number fluctuations of amplitude close to 1 should give rise to a GW 
signal as given by eq. ( |43| ) which may be detectable by LISA. Although this model does not provide 
an existence proof of GWs produced by neutrino inhomogeneous diffusion it gives, in our opinion, at 
least a reasonable plausibility argument. 

5.1 MHD effects 

In Ref. H it was showed that magnetic fields should be produced during the fluid acceleration process. 
This is a consequence of parity violation in the standard model which turns into a difference between 
the v a e~ and v a e + cross sections. It follows that an electric current J ex t appears in the e + e~ plasma 
in the presence of a net flux of neutrinos. Although the magnetic field produced by this current is 
initially very small, the coherent motion of the fluid induced by the neutrino flow amplifies this seed 
exponentially with time until equipartition is reached between the magnetic and fluid kinetic energies. 
Numerical simulations || show that in the range of parameters which is interesting from the point of 
view of GW detection, equipartition is reached well before fluid motion is damped by the viscosity 
and the universe expansion. It is clear that in such a situation the effects of the magnetic field on 
turbulence development have to be taken into proper account. 

The results we derived in Sec. ^ find here a natural application. Let us first observe that, since 
equipartition is reached already during the stirring phase, the equality us — va holds in this case. 
As we discussed in Sec. ||, depending on the characteristic of the energy injection spectrum P ex t(k) 
two cases have to be distinguished. If P e xt(&) is a decreasing function of k (7 < 2/3), or it is a delta- 
function peaked at k = k$ (7 = 0, as for Kolmogorov type turbulence), the turbulent velocity will be 
always smaller than the Alfven velocity. This means that MHD corrections can never be disregarded 
and the second of eqs. (pC| ) (k > A; oq ) applies for the turbulent velocity spectrum. If, however, 
-fext(^) is growing with k (7 > 2/3), or it is a constant (7 = 2/3, this is the case for a saturated 
spectrum), we always have > va and MHD effects are less important. 

6 GWs from turbulence produced by a first-order phase transition 

A first order phase transition proceeds through the nucleation of bubbles of the true-vacuum phase 
inside the false vacuum. Once a bubble is nucleated, if its radius is larger than a critical value, it 
begins to expand at velocity Vb, which for very strong phase transitions can approach the speed of 
light Once the bubble walls begin to collide, they break spherical symmetry and thus stir up the 

10 Note that fc cq = k$ in this case. 
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plasma at a scale comparable with their radii at the collision time. If the Reynolds number of the 
plasma is high enough, the energy released in coherent motions of the fluid is transferred to smaller 
scales and a turbulent spectrum establishes. It must be stressed that the stirring process must last 
for an enough long time in order to give rise to a fully developed turbulence. If this is not the case, as 
long as we are interested in GW production we can treat the resulting partially developed turbulence 
as if it were a fully developed one persisting for a time interval of the order of the turnover time on 
the largest stirring scale Q. 

The rate of nucleation of one critical bubble per unit volume, T, is suppressed by the exponential 
of the bubble Euclidean action. As the temperature of the Universe decreases the rate T gets larger 
and larger; when at time t* = the rate becomes comparable to H A the transition begins. In a 
neighbourhood of that time one can expand the Euclidean action, thus getting |i 

T-Hte^; (58) 

(3~ l sets the time scale of the process: in a time interval of order the transition is complete, 
i.e. all the volume has been converted to the true vacuum phase. The characteristic length scale of 
the process is thus Vb j3 . As pointed out in ref. [18], in realistic cases [5 is much larger than H*, so 
that cosmic expansion can be neglected for the whole duration of the process. 

We now want to understand what turbulent spectrum one can expect from such a first order phase 
transition. At generic time t the distribution of bubbles with a given radius R is |18| 



^(t) = -T(t R )p(t R ), (59) 
an Vb 

where p{t) is the probability that at time t a random point is in the false vacuum state and t R = t—R/vb 
is the time of nucleation of a bubble that at time t has radius R (the radius at time of nucleation 
is completely negligible). The factor p(tn) suppresses the bubble distribution at small radii, because 
smaller bubbles were nucleated at later times, when the fraction of volume still in false vacuum phase 
was smaller. We will show that collisions of small bubbles give a negligible contribution to the stirring 
of the plasma, so that from now on we will conservatively set p(tn) = 1. 

The energy carried by the expansion of a bubble of radius R is proportional to its volume R 3 . The 
energy distribution of bubbles at time t is thus ^(t) oc e' 34 R 3 e~ /3R / v *>, which in momentum space 
reads 

rl~F, 27rfl 

=-(t) oc f& k- 5 , (60) 

dk 

where we set k = 2tt/R. The plasma stirring process takes place when two (or more) bubbles collide: 
we expect that a collision will release energy in bulk motions at a scale comparable to the radius of 
the smaller colliding bubble, and the amount of energy released will be approximately proportional to 



17 



its energy. Since the energy distribution of small bubbles at any time t decreases like k~ 5 , from eq. (||) 
we expect to have no substantial deviation at large k from Kolmogorov spectrum when turbulence 
develops. On the other hand, in principle we can expect deviation from Kolmogorov spectrum at 
scales k ~ i, e , a t the characteristic scale of the transition. In order to understand if this is the 
case we must estimate the width of the stirring spectrum P ext (k): it is reasonable to assume that to 
be of the same order of magnitude of the width of the bubble energy distribution (|60|). The latter is 
peaked at 

k s = - — (61) 
5v b 

and has a width Ak ~ 1.4 x kg, if we define the width as the range of momenta k out of which the 
energy distribution dE/dk gets smaller than 1/e times its maximum value dE/dk\k s - 

The width Ak is small, i.e. most of the energy of the transition is carried by bubbles of characteristic 
scale very close to ks- This means that the stirring spectrum P ex t(k) has a width which is below (or 
comparable to) the precision of the results on turbulent motions, which are usually derived under the 
uncertainty condition Ak/k ~ 1 (see section ^). Therefore, we do not expect any substantial deviation 
from Kolmogorov spectrum in the case of turbulence produced by a first order phase transition, i.e. we 
can set 7 = in all the formulae of section |I[ 

It is interesting to deal with the electroweak phase transition (EWPT), because it takes place at 
the temperature ~ 100 GeV and therefore, if it is of the first order, it gives rise to a GW background 
peaked at a the frequency (see eqs. ( fl4| ) and (61)) 
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(62) 



which is in the range of maximum sensitivity of LISA if us is relativistic and (3/H* ~ 100. As pointed 
out in ref. the duration of the turbulent source of GWs is roughly the turnover time on the largest 
length scale, At* ~ P~ 1 Vb/us- The amount of relic GWs at milliHertz frequencies is thus (see eq. (|45|) ) 
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(63) 



If the bubbles expand as detonation fronts, i.e. faster than sound, the expansion velocity of the 
bubbles, Vb, and the turbulent velocity on the largest scale, us, can be derived in terms of a, the ratio 
between the false vacuum energy and the plasma thermal energy at the transition time || B]. One 
obtains 

(a 2 + 2a/3) 1/2 



v b (a) 
u s {a) 



1 + a 



a 



0.72 a + 4/27 
1 + 0.72 a 
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(64) 
(65) 



plugging these velocities into eq. ( 62 ) and taking the weak transition-small a limit one gets the result 
of ref. 0, apart from the factor of order 4 discussed in section f|. This turns out to be a signal too low 
for LISA. Clearly a more interesting signal can be reached in the strong transition-large a limit: this 
case could be motivated by supersymmetric extensions of the Standard Model, which in some small 
regions of parameter space give rise to a very strong first order EWPT. This strong transition, apart 
from creating a large amount of relic GWs, can be responsible for the generation at the electroweak 



scale of the observed baryon asymmetry in the Universe (see Ref. [19] for a review). An extensive 
numerical study of such EWPTs in the Minimal Supersymmetric Standard Model and in the Next-to- 
MSSM has been performed in ref. ||. If the transition is very strong a is close to 1 and both velocities 
Vb and us reaches relativistic values; the intensity of relic GWs from turbulence thus becomes 
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The stronger is the transition and the smaller is (3/H*, because the system experiences a larger degree 
of supercooling. However, for a strong EWPT (3/H* ~ 100 is not unusual (see ref. |l^] for analytic 
estimates and ref. p] for exact numerical computations in specific examples). Such values of (3/H* 
would give a signal above the planned sensitivity of LISA. 



6.1 MHD effects 

So far we ignored possible MHD effects which may arise because of the presence of magnetic fields 
during the EWPT. This approach, however, is not always justified since magnetic fields may have 
been produced before or during the EWPT by a number of different processes (for a review see [jlQ]). 
Quite solid arguments in favor of magnetic field generation at a first order EWPT were presented by 
Baym, Bodeker and McLerran ||. Magnetic fields are seeded by the motion of the dipole charge layer 
which is produced across the bubble boundaries because of the small excess of top quarks over top 
anti-quarks and of the potential barrier which quarks face at the bubble wall. Even if the seed fields 
are initially very tiny, they are amplified exponentially with time by the strong turbulence which is 
produced when bubble walls collide. This process is very effective due to the high bubble wall velocities 
and the high value of the Reynolds number at that time, Re ~ 10 12 ||. The amplification processes 
ends when energy equipartition is reached between the magnetic field and the turbulent velocity field. 
We will show now that the result of Baym et al. may have relevant consequences for the GW signal 
to be expected at LISA. 

We apply here the results we derived in sec. ||. We have shown there that the MHD turbulent 
energy spectrum takes the same form as for the Navier-Stokes type of turbulence upon a simple 
rescaling of the parameter 7 (defined in eq. ([?])) and of the wavenumber k$ (see Eqs. (p2|)). Due to 
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Figure 1: The figure shows whether an EWPT characterized by a pair (a, (3) gives rise to a GW back- 
ground from turbulence which is detectable by LISA. The shaded region is "detectable"; the dashed 
and continuous lines refer respectively to the Navier-Stokes and to the Magneto-Hydro-Dynamical 
case. 



the huge Reynolds number we think it is a reasonable approximation to assume that equipartion is 
reached already at the maximal stirring scale, i.e. k eci = ks- Furthermore, since we showed above that 
7 = for the EWPT with no magnetic fields, in the MHD we have just to replace 7 with 5 = 1/6. By 
replacing this value in eq. (]4l|) and following the same procedure we used to derive eq. (|63|), we get 
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(67) 



This is a quite interesting result! For instance, by taking a ~ 1 the velocities (p3) and ( pop becomes 
relativistic and the overall numerical factor in eq. ( |6"7| ) becomes 1.8 x 10~ 12 : by comparing it with 
eq. (|66|) we see that, even taking into account the smaller enhancement due to the 2/3 exponent for 
(3/H*, MHD effects amplify the signal to be expected at LISA by one order of magnitude. This signal 
is certainly within the planned sensitivity of the instrument Q. The strong enhancement with respect 
to eq. ([3f]) is easily explained by the milder slope of the MHD turbulent energy spectrum at high 
frequencies with respect to the Kolmogorov's. 

By fixing the temperature at which EWPT happens at 100 GeV, eqs. (63) and ( p7| ) are functions 
of a, (3 and /; the LISA (expected) sensitivity in ti^O,^ clearly depends only on frequency / and it 
is roughly 10 -11 at 0.1 mHz, 10~ 12 at 1 mHz and 10~ n at 10 mHz (see Q for the LISA sensitivity 
curve). Figure [i] shows the values of a and (3 for which there exists at least a range of frequency in 
which the signals (|63| ) and ( |67|) are above the sensitivity of LISA. Clearly large values of a, i.e. strong 
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transitions, are favoured. Notice that too large values of (3/H* are disfavoured, because in such cases 
the characteristic frequency ( |62|) is too high for LISA. Notice the effect of the MHD corrections, which 
substantially enhance the detectability of the signal mostly at small values of f3/H*, because in such 
cases the characteristic frequency fl62| ) is quite below the sensitivity window of LISA. 

Furthermore, we should keep in mind that tangled magnetic fields act themselves as a source of 



GWs [13]. Assuming that magnetic fields are in equipartition with the turbulent velocity field, this 
source should enhance the expected GW signal by a factor of two in all the frequency range over which 
equipartition is established (see eq. (|2l|)). 



7 Conclusions 

In this paper we investigated the characteristics of the GW signal sourced by turbulence in the early 
universe. We extended previous work on the subject by considering the effects on the turbulence 
energy spectrum of a continuous stirring power and of magnetic fields. In the first part of our work 
we derived a suitable parameterization for the turbulent energy spectrum in the case of a continuous 
stirring power with and without external magnetic fields. Interestingly we found that the MHD tur- 
bulent spectrum can be written in the same form as the Navier-Stokes ones upon a simple rescaling 
of only two parameters. Our treatment is original and may find applications which go beyond the 
aims of this paper. Following the approach of Kosowsky et al. Q, who consider only Kolmogorov type 
of turbulence, we used our formulation to determine the characteristics of the GW signal produced 
by cosmological turbulence with a more general energy spectrum. In the Kolmogorov case we found 
minor discrepancies with the result reported in 0]. Deviations from the Kolmogorov spectrum may 
have relevant consequences for the GW signal. Indeed, since the modified spectra are generally less 
steep than the Kolmogorov's, these deviations will result in a stronger signal at high frequencies. Fur- 
thermore, our analysis may allow to extract valuable informations about the nature of the turbulence 
source and the presence of primordial magnetic fields from the GW background power spectrum which 
may be measured by forthcoming experiments. 

We applied our results to estimate the GW expected signal for two possible kind of sources. In 
sec. H we considered GW production by neutrino inhomogeneous diffusion according to the mechanism 
proposed in Ref. ||. We showed that a detectable signal can be produced only if the amplitude of 
the lepton number fluctuations is close to unity over a wide range of wavenumbers. This possibility 
is not unreasonable as active-sterile neutrino oscillations or some models of leptogenesis based on the 
Affleck-Dine scenario can indeed give rise to domains with opposite lepton number. Although we 
argued that the former scenarios can hardly give rise to GW signal detectable by LISA, the latter 
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scenario offers more promising observational perspectives. 

Finally in sec. ^ we considered GW production by turbulence at the end of a first order phase 
transition. In the absence of strong magnetic fields no substantial deviations have to be expected from 
the results of Ref. H based on the assumption of a Kolmogorov turbulence spectrum. In the case of 
the EWPT the GW signal may be above the LISA planned sensitivity only if the transitions is very 
strong with large bubble wall velocities. This may be achieved in next to minimal extensions of the 
supersymmetric standard model ||. A more favorable situations may be obtained if strong magnetic 
fields were present at the end of the transition. This is actually to be expected in the case of the 
EWPT if it is first order We showed that in this case the signal can be strongly enhanced with 
respect to the non-magnetic case and be detectable by LISA if bubble wall velocity was not too much 
smaller than unity. Similar consideration may apply to other physical situations like, for example, 
during reheating at the end of inflation. These results open new perspectives for a successful detection 
of GWs backgrounds produced in the early universe. 
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Appendix 

Computation of the spectrum /ig^gw(/) 

In this appendix we work out the detailed computation of the relic spectrum /iqH (J) in the three 
cases 7 < 2/3, 7 ~ 2/3 and 7 > 2/3, along the strategy described in section |4j. 
If 7 < 2/3 eq. (H) becomes 

( Mk, *.) W Q > * 27 'p^$ W)2 At, 4 fcf" 7 * 3 (k " k') , (68) 

which in real space reads 

( h^x, Q %(x, Q ) ~ 27- V2 tt 2 (Gwf ^ ^ ^ j k » fc -i4/3+ 7dfe (6g) 

2 — 37 J kg 

- ^^^v^r^y 1 ^?. (70) 



2-37 " " A, \k s 



Changing variable from the turbulent scale k to the characteristic frequency / (see eq. (34)), we have 
( MM*) Mx,t.) ) - ( 2 -3 7 )(4 + 37) ^* Us k S J fs {j-J J ' 
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where fs = fk s is the characteristic frequency of the largest length scale. The characteristic amplitude 
squared is thus 
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The characteristic amplitude squared measured today at a frequency / (see eq. (pS|) ) is 
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where we used Gw = H*/(2ir) and both frequencies / and /g have been redshifted. By means of 
eqs. ( ^Q) and (f44[) we finally have the spectrum, eq. (41). 

The above computations can be easily generalized to the 7 > 2/3 case: 
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The characteristic amplitude measured today is 
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which, combined with eq. (|40|), gives eq. 

Looking at eq. (|3l|) we see that the intermediate case 7 ~ 2/3 is formally the same as 7 < 2/3 case, 
with the replacements 7 — > 2/3 and 2 _3 7 — > 1, and with an extra logarithmic factor. The present-day 
observable are thus 

ksr 3 \ n r 2 r 9* r 2/ V/ x ~ 3 
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/s 
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(78) 



and the spectrum given in eq. ([43]) . 

As a cross-check, notice that the 7 ~ 2/3 case is also equal to the 7 < 2/3 one, again modulo 
proper replacements and the logarithmic factor. The present-day observable computed in the two 
ways are equal. 
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